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We formulate a kinetic theory for non-centrosymmetric superconductors at low temperatures 
in the clean limit. The transport equations are solved quite generally in spin- and particle-hole 
(Nambu) space by performing first a transformation into the band basis and second a Bogoliubov 
transformation to the quasiparticle-quasihole phase space. Our result is a particle-hole-symmetric, 
gauge-invariant and charge conserving description, which is valid in the whole quasiclassical regime 
■ (lq| <C kp and hcu <^ Ep). We calculate the current response, the specific heat capacity, and 

OA ' the Raman response function. For the Raman case, we investigate within this framework the 

, polarization-dependence of the electronic (pair-breaking) Raman response for the recently discov- 

er ■ ered non-centrosymmetric superconductors at zero temperature. Possible applications include the 

systems CePtaSi and Li2Pdi:Pt3_a;B, which refiect the two important classes of the involved spin- 
orbit coupling. We provide analytical expressions for the Raman vertices for these two classes and 
calculate the polarization-dependence of the electronic spectra. We predict a two-peak structure 
and different power laws with respect to the unknown relative magnitude of the singlet and triplet 
contributions to the superconducting order parameter, revealing a large variety of characteristic 
fingerprints of the underlying condensate. 



o 



o 



X 



INTRODUCTION 



Oh: 

■ In a large class of conventional and in particular unconventional superconductors a classification of the order 
J parameter with respect to spin singlet /even parity and spin triplet/odd parity is possible, using the Pauli exclusion 
principle. A necessary prerequisite for such a classification is, however, the existence of an inversion center. Something 
of a stir has been caused by the discovery of the bulk superconductor CePtaSi without inversion symmetry 
which initiated extensive theoretical 12 , 2g and experimental studies 0, [llj . In such systems the existence of an 



C antisymmetric potential gradient causes a parity-breaking antisymmetric spin-orbit coupling (ASOC), that gives rise 
^ to the possibility of having admixtures of spin-singlet and spin-triplet pairing states. Such parity-violated, non- 
I— I, centrosymmetric superconductors (NCS) are the topic of this chapter, which is dedicated particularly to a theoretical 
study of response and transport properties at low temperatures. We will use the framework of a kinetic theory 
I described by a set of generalized Boltzmann equations, successfully used before in to derive various response and 
J — • transport functions such as the normal and superfiuid density, the specific heat capacity (i. e. normal fraction and 
, condensate properties, that are native close to the long wavelength, stationary limit) and in particular the electronic 
Tj" I Raman response in NCS (which involves frequencies hv comparable to the energy gap Ak of the superconductor). 

[ A few general remarks about the connection between response and transport phenomena are appropriate at this 
CO ' stage. Traditionally, the notion of transport implies that the theoretical description takes into account the effects of 
[ quasiparticle scattering processes, represented, say, by a scattering rate F. Therefore, we would like to demonstrate 
with a simple example, how response and transport are intimately connected: consider the density response of normal 
metal electrons to the presence of the two electromagnetic potentials and A'^'*^*, which generate the gauge- 

invariant form of the electric field E = — V<I>'"^* — dA/cdt. In Fourier space (V — iq, d/dt — > —iuj) one may write for 
the response of the charge density: 



' Sue = e^iq • Mo(q, to) ■ E 

with Mn the Lindhard tensor and q • Mq • q = Mq the Lindhard function, appropriately renormalized by collision 
effects ii: 



Mo(q,w) - 



1 - 



1 _ Ca{(i,uj+ir) 
£o(q,0) 



uj+ir 

Here £o('?,w) denotes the unrenormalized Lindhard function in the collisionless limit F — > 0: 



-i-oiqjWj — TT > r— 

V ^ ep+q/2 - ep-q/2 " nuj 
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In this definition of tlie Lindliard function, denotes the equilibrium Fermi-Dirac distribution function and = 
Ck+M represents the band structure. Now the aspect of transport comes into play by the observation that Mo(q, uj+iT) 
may be expressed through the full dynamic conductivity tensor cr(q, w) = e^(9n/9/^)D(q, w) of the electron system 
as follows: 

^o(q,w = — ^ — : 

lU) ~ q - 13(q, u!) • q/(l — iujt) 

with q • (T • q = iwe^£o(q, i^) and with the so-called diffusion pole in the denominator of Mo(q, cj) reflecting the 
charge conservation law. This expression for the Lindhard response function Mq clearly demonstrates the connection 
between response (represented by Mq itself) and transport (represented by the conductivity cr), which can be evaluated 
both in the clean limit F — !■ and in the presence of collisions F 7^ 0. In this sense, the notions of response and 
transport are closely connected and therefore equitable. In this whole chapter we shall limit or considerations to the 
coUisionless case. 

An important example for a response phenomenon involving finite frequencies is the electronic Raman effect. Of 
particular interest is the so-called pair-breaking Raman effect, in which an incoming photon breaks a Cooper pair 
of energy 2Ak on the Fermi surface, and a scattered photon leaves the sample with a frequency reduced by 2Ak/^, 
has turned out to be a very effective tool to study unconventional superconductors with gap nodes. This is because 
various choices of the photon polarization with respect to the location of the nodes on the Fermi surface allow one 
to draw conclusions about the node topology and hence the pairing symmetry. An example for the success of such 
an analysis is the important work by Devereaux et al. 6] in which the d2;2_y2-symmetry of the order parameter in 
cuprate superconductors could be traced back to the frequency-dependence of the electronic Raman spectra, that 
directly measured the pair-breaking effect. Various theoretical studies of NCS have revealed a very rich and complex 
node structure in parity-mixed order parameters, which can give rise to qualitatively very different shapes, i. e. 
frequency dependencies, of the Raman intensities, ranging from threshold- and cusp- to singularity-like behavior. 
Therefore the study of the polarization dependence of Raman spectra enables one to draw conclusions about the 
internal structure of the parity-mixed gap parameter in a given NCS. 

This chapter is organized as follows: In section □ we introduce our model for the ASOC, the two order parameters 
on the spin-orbit splitted bands and the pairing interaction. Then, in section □ we derive the kinetic transport 
equations for NCS at low temperatures in the clean limit and transform these equations into the more convenient 
band-basis. In section „ the transport equations are solved quite generally in band- and particle-hole (Nambu) space 
by first performing a Bogoliubov transformation to the quasiparticle-quasihole phase space and second performing the 
inverse Bogoliubov transformation to recover the original distribution functions. We demonstrate gauge invariance 
of our theory in section □ by taking the fluctuations of the order parameter into account. Within this framework, we 
calculate the normal and superfluid density in section □ and the specific heat capacity in section „ In section „ our 
particular interest is focussed on the electronic Raman response. We investigate the polarization-dependence of the 
pair-breaking Raman response at zero temperature for two important classes of the involved spin-orbit coupling. 
Finally, in section □ we summarize our results and draw our conclusions. 



ANTISYMMETRIC SPIN-ORBIT COUPLING 



We start from a model Hamiltonian for noninteracting electrons in a non-centrosymmetric crystal [25j 

H = ^cl^ i^-kSaa' + 7k ■ "^fTtr'] Cka' , (1) 

where ^k represents the bare band dispersion assuming time reversal symmetry (^_k = Ck), cr, cr' —'[,■1 label the spin 
state, and t are the Pauli matrices. The second term describes an antisymmetric spin-orbit coupling (ASOC) with a 
(vectorial) coupling constant 7k. The pseudovector function -y^ has the following symmetry properties: 7_k = ^7k 
and 57g-ik = 7k- Here g denotes any symmetry operation of the point group G of the crystal under consideration. 
In NCSs two important classes of ASOCs are realized, reflecting the underlying point group G of the crystal. In 
particular, we shall be interested in the tetragonal point group C4,v (applicable to the heavy Fermion compound 
CePtsSi with Tc=0.75 K Q for example) and the cubic point group O (applicable to th e sy stem Li2Pd2;Pt3_a;B with 
Tc=2. 2-2.8 K for x=0 and Tc=7.2-8 K for x=3 Q). For G = C^y the ASOC reads 



7k = g±{k X e^) + g\\k^kyk^{kl ~ kl)e^ . (2) 



(a) 



(b) 
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FIG. 1: The angular dependence of |7^| for the point groups C41, and O. Since dk||7it, these plots show also the magnitude of 
the gap function in the pure triplet case for both point groups. 



In the purely two-dimensional case (g|| — 0) one recovers, what is known as the Rashba interaction 0,i[ii. We will 
choose for simplicity g|| ~ for our Raman results. For the cubic point group Q = 7^ reads (30l | 



7k = gik - g3 k^{kl + fc2)e^ + ky{kl + kl)&y + k:,(kl + /c^)e^ 



(3) 



where the ratio gs/gi ~ 3/2 is estimated by Ref. [30]. Because of the larger prefactor g3 > gi we will keep the higher 
order term for our further considerations. Thus, in terms of spherical angles, k = (cos(/)sin0, sin sin 0, cos0), the 
absolute value of the 7j(.-vectors for both point groups, illustrated in Fig. [1] reads 



l7kl = sin 61 



l7kl = \/l-^sin^2^-Asin4, 



sin^ 20 (9 sin^ ( 



-4 



for 6*4^ 
for O 



(4) 

(5) 



By diagonalizing the Hamiltonian [Eq. ([T])], one finds the eigenvalues ^A(k) = ■^k + -^iTkL which physically corre- 
sponds to the lifting of the Kramers degeneracy between the two spin states at a given k in the presence of ASOC. The 
basis in which the band is diagonal can be referred to as the hand basis where the Fermi surface defined by C±(k) = 
is splitted into two pieces labeled ±. Sigrist and co-workers have shown that the presence of the ASOC generally 
allows for an admixture of a spin-triplet component to the otherwise spin-singlet pairing gap 12]. This implies that 
we may write down the following ansatz for the energy gap matrix in spin space: 



,(k) = {[^k(r)l + dk(T) 



(6) 



where V'k(T') and dk(T) reflect the singlet and triplet part of the pair potential, respectively. In the band basis we 
find immediately 



A±(k) = 7/;k(T)±|dk(T)| 



(7) 



It has been demonstrated that a large ASOC compared to k^T^ is not destructive for triplet pairing if one assumes 
dk||7k ESi: 



dk(T) = d{T)% , (8) 

whereas the temperature-dependent magnitudes i^iT) and d(T) of the spin-singlet and triplet energy gaps are solu- 
tions of coupled self-consistency equations and 7k is defined by 

7k 



7k = 



(9) 
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Thus the energy gap of Eq. (O can be written as: 

A±(k)=V(T)±d(r)|7k| . (10) 

For the T = Raman response in section □ we will use the following ansatz for the gap function on both bands 
(+ and -) [13: 

A±(k) = 4,± d|7k| = (1 ± p|7kl) = A± , (11) 

where the parameter p — d/^p represents the unknown triplet-singlet ratio. Accordingly, the Bogoliubov-quasiparticle 
dispersion is given by i?±(k) ~ ^±(k) + Aj_(k). If we assume no q-dependence of the order parameter, AA(k) [and 
also _E;!^(k)] is of even parity i.e. A;^(— k) = AA(k). It is quite remarkable that although the spin representation of the 
order parameter Ao.CT'(k) has no well-defined parity w.r.t. k — ^ — k, as easily seen in Eq. (jS]), the energy gap in band 
representation has. Note that for Li2Pd2,Pt3_j.B the parameter p seems to be directly related to the substitution of 
platinum by pall adium, since the larger spin-orbit coupling of the heavier platinum is expected to enhance the triplet 
contribution VM- This seems to be confirmed by penetration depth experiments (sol . Isij . 
The corresponding weak-coupling gap equation reads 

A,(k,T) = -5]y4^A^(k',r)0^(k') (12) 

k',/j. 



with 



and its solution are extensively discussed in Ref. [ISf . Here and in the following we choose a separable ansatz for the 
pairing-interaction (cf. Ref. with Cm = 0, i.e. without Dzyaloshinskii-Moriya interaction): 

^kk' = r, + A/ir,|7kll7k'l , (14) 

where T s and Tt represent the singlet and triplet contribution, respectively. Although an exact numerical solution 
of Eqs. p2|) - (IT4l) with a microscopic pairing interaction would be desirable, we restrict ourselves in this work to a 
phenomenological description which allows an analytical treatment of response and transport in NCS. 



DERIVATION OF THE TRANSPORT EQUATIONS 

In this section, we study the linear response of the superconducting system to an effective external perturbation 
potential of the form 

S^Taa' = [e$(q,^) - • A(q,cj)] 5,,, + £^Af(q,cj)-^^2^k^Af (q, w)5,,, . (15) 

Here $ and A denote the electromagnetic scalar and vector potential. Electronic Raman scattering is described 
in addition by the third term in Eq. (1151) . It describes a Raman process where an incoming photon with vector 
potential A^, polarization and frequency a;/ is scattered off an electronic excitation. The scattered photon with 
vector potential A'^, polarization e'^ and frequency u)s = ijJi — uj gives rise to a Raman signal (Stokes process) and 
creates an electronic excitation with momentum transfer q. Further, the Raman vertex in the so-called effective-mass 
approximation reads 

In general, an external perturbation can be decomposed into a vertex function Ok and a related potential 5£^a ■ 

'^Cr=E°k'5?a. (17) 
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TABLE I: External perturbations can be decomposed into a vertex-function and a potential. The vertex-function is charac- 
teristic for each response-function and can be classified according to parity (w.r.t. k — >■ — k) and dimension. 



vPTtp"y 
Ok 


(Active) potential parity 


d 1 vn PTi nti 




e 

evk 

( Ex(k) ei5A(k)\ 
\ T dT J 


^cxt 
^ext 

ST 


even 

odd 

even 


scalar 
vector 
scalar 


charge and 
current 
specific heat 
capacity 


m(M-i),,, 


roAfAf 


even 


tensor 


Raman 



A list of all relevant vertex functions and potentials that will be discussed in this chapter is given in Table HI The 
charge density response to the electric field E = — V^^'^* — dA/cdt is characterized by a constant vertex at = e 
(electron charge) and therefore of even parity (w.r.t. k — — k), whereas the current response to the vector potential 
A depends on the odd vertex-function ak = eVk (electron velocity). In case of the specific heat capacity Cy(r), the 
role of the fictive potential is played by the temperature change 5T, which couples to the energy variable ^k- For the 
Raman response, this fictive potential depends essentially on the vector potential of the incoming and scattered light. 
The response and transport functions will be obtained as moments of the momentum distribution functions with the 
corresponding vertex (see section n-J. 

In addition to the external perturbation potentials, also molecular potentials can be taken into account within a 
mean-field approximation: 

= <5er + E (/kp + v^) Sfip = "k-^Ca + v^sm + J2 /kp'^np . (18) 

per a pa 

The short-range Fermi-liquid interaction f^^ leads to a renormalization of the electron mass [IJ] and the long-ranged 
Coulomb interaction with Vq = Ane^/q^ is included self-consistently through the macroscopic density fluctuations 
6ni = J2pa '^"'P with the non-equilibrium momentum distribution function Sup. 

The potentials S£,^* are assumed to vary in time and space esc exp{iq^-r—iujt). Then the response to the perturbation 
potentials can generally be described by a nonequilibrium momentum distribution function n^p, , which is a matrix in 
Nambu, momentum and spin space with p = k + q/2 and p' = k — q/2. The evolution of the nonequilibrium matrix 
distribution function in time and space is governed by the matrix-kinetic (von Neumann) equation 0, [29| 

^^ipp' + E[^pp"'ip"p']_ = ° (^9^ 
p" 

in which the full quasiparticle energy plays the role of the Hamiltonian of the system. This equation holds for 
hio ^ Ef and |q| <^ kp. In general, a collision integral (see e.g. [2^) could be inserted on the right hand side of 
Eq. ([T^ that accounts for the relaxation of the system into local equilibrium through collisions. In the following we 
will assume the absence of collisions 32] . After linearization according to 

nip"p' = 2ik(q.'^) =2ik'^q,o + '5nk(q,w) (20) 

ip"p' = ik(q' ^) = ik^i-o + '^ik(q' ^) ' (21) 

the matrix-kinetic equation assumes the following form in cj — q- and spin-space: 

h^Sn^ + 6n^^l_ - ^l^Sn^ = S^^nL - n^Ji^ . (22) 

Here, uj is the frequency and k± = k ± q/2, with q representing the wave number of the external perturbation. The 
equilibrium distribution function and quasiparticle energy are matrices in Nambu and spin space: 



nk gk 

gi 1- n_k 



(23) 



K -[^k + 7k 
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The momentum and frequency-dependent deviations from equilibrium are defined as 

Sng 

and 



(25) 



(26) 



respectively. In the spin basis, the matrix-kinetic equation [Eq. (|22p] represents a set of 16 equations, which can be 
reduced to a set of 8 equations by an unitary transformation into the band basis (also referred to as helicity basis). 
This SU(2) rotation is given by [28| 



Uk 

u* 



exp 1 -*yn^ 



= cos ■ 



m-, • T sm — - 

^ 2 



7k X z 
l7k X z| 



(27) 
(28) 
(29) 



which corresponds to a rotation in spin space into the z-direction about the polar angle 6-y between 'yk ^^'^ ^- Then 
Eq. (p2|) may be written as 



^IZk+'^ZikiZk- 



iki^k- 



= t[U'5e^c/k-C£L2ik-iik- -ciUHk+tik+cilc+^ekC^k- 



or, more simply 



where the equilibrium distribution function and energy shifts in the band basis are given by 






V -a;6 












i(l + e-^-) 








(30) 



(31) 



(32) 



i(i + e+e+)y 



and 



r = 






Va; 







-Al 






-A_ 




A+ 





(33) 



The deviations from equilibrium can be parameterized as follows: 



^ik 



I£[+'^2ikIZk- 



/ 5n\_ 








\ 





S-nt 










-Sc/_* 


-Snt 













-Sn\_ J 


( 








SAX \ 







-SAt 








-(JA^L* 







\ 5A^* 











(34) 



(35) 



Thus, we have now derived a set of equations in spin- and band-basis [Eqs. p2|) and ((3T|) ] that allow us to determine 
the diagonal and off-diagonal non-equilibrium momentum distribution functions. In sections □-□ we will use these 
distribution functions to determine the normal and superfluid density, specific heat capacity, and the Raman response 
of NCS. From now on, we will omit the index "&" indicating the band-basis, since all further considerations will be 
made in the band-picture. 
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SOLUTION BY BOGOLIUBOV TRANSFORMATION 



In what follows we will solve the kinetic equation pip , derived in the previous section. For this purpose, we per- 
form first a Bogoliubov transformation into quasiparticle space, where the kinetic equations are easily decoupled and 
then solved. For the subsequent inverse Bogoliubov transformation we will introduce parity projected quantities to 
obtain finally a relation between the diagonal and ofi^-diagonal energy-shifts on the one side and the non-equilibrium 
distribution functions on the other. As a fist step towards the solution of the kinetic equations, the momentum distri- 
bution matrix and the energy matrix (both in band-basis) are diagonalized through the following Bogoliubov 
transformation 



f f{E+) \ 

f{E^) 

f{-E-) 

V fi-E+) J 

f E+ \ 



-E- 

\ -E+ J 



(36) 



(37) 



with the Fermi-Dirac distribution function f{E\) — [exp{E\/kBT) + 1] ^. The Bogoliubov matrix has been found to 
read in the band basis 



/ u+ w+ \ 

u- -V- 

v*_ u- 

\ ~vX u+ J 



(38) 



with the coherence factors 



'1 

2 \ Exik) 



/I A gA(k)\ AA(k) 

E,{k)J\A,{k)\ 



(39) 
(40) 



satisfying the condition lii^l 



1, by which the fermionic character of the Bogoliubov quasiparticles is 



established. In order to solve the transport equation in the band basis (|3ip. one may multiply from the left with the 
Bogoliubov matrix B_l^^ and from the right with . The result is 



or, more simply 



kiik- 



- St+'^ik^k-^j^-lik-^k- -^U«k+^k+s],+<5ek^k- 



(41) 



(42) 



The new Bogoliubov-transformed quantities describing the deviation from equilibrium are identified from the preced- 
ing equations and labeled as follows: 



/ dv+(k) (57+ (k) \ 

(5i^_(k) -(57_(k) 

-(57l(k) -(5i^_(-k) 

~Siy+{-k) J 

/ 6E+{k) SD+{k) \ 

SE^{k) -SD^{k) 

-SD*_{k) -SE^{-k) 

\(5L>;(k) -6E+{-k) I 



(43) 



(44) 
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The solution of Eq. (|42)) for the quasiparticle distribution functions is the set of the following eight equations (A = ±) 



<5i/A(k) 
<57A(k) 



^~Vx (k) 



uj + Vx (k) 



yx{k)6Ex{~k) 



u;-77+(k) 
^;^(k) 



exik)SDl{k) , 



W + 7/+(k) 

where we have introduced the following abbreviations: 

ritik) = Exik+)±Exik^) , 



and 



yA(k) 



eA(k) = 



_f[Ex(k+)]^nEx{k-^ 
Ex{k+) - Ex{k-) 

1 - nEx{k+)] - f[Ex{k-)] 
Exik+) + Ex{k-) 



(45a) 
(45b) 
(45c) 
(45d) 

(46) 
(47) 

(48) 



The expressions for these quantities in the long-wavelength limit can be found in appendix 1. In this limit, the 
difference quotient yx(k) is equal to the Yosida kernel yA(k) which is given by the derivative of the quasiparticle 
distribution function 



yA(k) 



df[Ex{k)] 
dExik) 



1 



1 



(49) 



and is crucial for the temperature dependence of all response and transport functions. Accordingly, 8a (k) '^^^ Sx{k) 
represents the kernel of the self-consistency equation ([T^ . It is instructive to note that the distribution functions 
Si^xO^) and (57a (k) have a clear physical meaning: The diagonal component Svxik) — 6{a''^ax){k) describes the 
response of the Bogoliubov quasiparticles (with the quasiparticle creation and annihilation operators a\, ax hi the 
band A). The off-diagonal component (57a (k) = S{axax){k) describes the pair-response. Note that the abbreviations 
Tj^ik) are of even (+) and odd (— ) parity w.r.t. k — k and become very simple expressions in the small wavelength 
limit (see appendix 1). 

For the inverse Bogoliubov transformation it is convenient to introduce parity-projected quantities which are labeled 
by s = ±1: 



Sn^'^k) = ^[Sn{k) + sSn{-k)] , 

S^<^^\k) = ^[Sm + sSa-k)] . 
In almost the same manner also the off-diagonal components are decomposed by 



and 



/^Hk) = i 



(5A(^)(k) = i 



SA(k) 



|A(k) 
A*(k) 



|A(k)| 
A(k) 



(5A(-k) 



(50) 
(51) 

(52) 

(53) 



|A(k)| |A(k)| 

We use the same symmetry classification for the Bogoliubov transformed quantities. The physical meaning of 
6Ax{k,q,uj) becomes clear after a decomposition into its real and imaginary part 

a(k, q,tj)e''^('i'") 

A(k) 



(5A(k, q,uj) 



A(k) 

[<5a(k, q,w)+i,5</.(q,w)|A(k) 



(54) 



|A(k)| 
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With Eq. (|53l) we can identify (5A(+^ (k, q, oj) = 5a(k, q, oj) as the amplitude fluctuations and ^ (k, q, a;)/A(k) = 
i6ip{q, uj) as the phase fluctuations of the order parameter. 

The off-diagonal energy shift i5A^''-'(k) can be determined from a straightforward variation of the self-consistency 
equation 



5A(^)(k) = ^y,\^<5g(^)(k') 

k'/j. 



(55) 



with Sg['\k) 



~6x{'k)5A''^\'k). This off-diagonal self-consistency equation will play an important role for the gauge 
invariance of the theory, as will be discussed in section „ 

From the symmetry-classification we can assign to each transport and response function (see Table|l| the correspond- 
ing momentum distribution function Sn^^\]<.) or Sn\^ ''(k): The vertex function of the (charge) density- and Raman- 
response is even in k, thus only the even distribution function 5n^^\'k) contributes to those response-functions. For 
the current-response (dynamic conductivity), the vertex-function (oo-(k) = evk) is odd in momentum. Thus, only 
(5n^ ^(k) contributes to the conductivity upon summation over k. Furthermore, the Bogoliubov-transformation can 
now be written in this simple form 



^7a 



(k) 



91^^ (k) P^^\k) 
-pi'\k) q^^\k) 



SE['\k) 



SDi^\k) 



9l^^(k) 



(k) 

-p['\k) q['\k) 



is) 

Px 



which might easily be inverted by using the sum rule 




^d^^(k) 
<5A('')(k) 



(56) 
(57) 



qi'Hk) + p\'>ik) 



= 1 



Here, we have defined the real- valued coherence-factors 

qi'\k) - \ux{k+)ux{k-)\ - s\vx{k+)vx{k~)\ 

and 

q['\k) = \ux{k+)vx{k~)\ + s\ux{k~)vx{k+)\ 

with the explicit form 



(58) 



(59) 



(60) 



and 




eA(k+)6(k-)-s|AA(k)|2 



2Ex{k+)Ex{k^ 



(61) 



I 2 2Exik+)Exik-) 
From Eqs. ((57| and we finally obtain the following solution of the matrix-kinetic equation 



(62) 



/ Sn+{k) \ 


1 


(5n^(k) 










V 



'11 



N2i 



Ni2 N: 

N21 N22 N: 

N31 N32 N33 
Nil N42 N, 



J 



( ^eA+(k) \ 

<5A+(k) 

V'^A,-(k)y 



(63) 



The vector on the left hand side contains the non-equilibrium momentum distribution functions [defined in Eq. (j34l) ] 
which can be expressed in terms of the diagonal and off-diagonal energy-shifts [defined in Eq. ((35)) and obtained from 
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Table H] and Eq. ([55]) ]. The matrix-elements Nij read in detail: 



Nil 

Nl2 
Ni3 



-pi+)^(k)ei+)(k) 



qi+^ {k)q\-> ik)y{-' (k) + p\^> {k)p{-> (kjO^^ (k) 
q^+\k)p[+\k) [yi+)(k)-el+)(ky 
qi+\k)pi-\k)y^-\k) q^-\k)p^+\k)ei-\k) 



(64a) 
(64b) 
(64c) 
(64d) 



N22 
N23 
N2A 



q^-^\k)y^^\k)+p^~^\k)Q^^\k) 
q'^-\k)p^^\k)y[-\k) - qt\k)p^-\k)Q^-\k) 
q[-\k)p^-\\.) [yi+)(k)-ei+)(k)" 



(64e) 
(64f) 
(64g) 



7V33 = pi+^'(k)yl+)(k) + ql+)2(k)ei+)(k) 

iV34 = pi+Hk)pi-)(k)yi-)(k) + gi+)(k)<zi-)(k)ei-)(k) 



(64h) 
(64i) 



7V44 - pi-^'(k)yW(k) + gr^'(k)er^(k) 



(-)2. 



(64j) 



The matrix-elements Nij are symmetric, i.e. Nij = Nji and the occurring products of coherence-factors can be found 
in the appendix 1. Above, we have defined the following abbreviations: 



yi^^(k) 



ei^^(k) 



Vx (k) 



^'-vi'^'ik) 

^-^i^^'(k) 



-yA(k) 



(65) 



-eA(k) 



The matrix-elements A^is, N23 and N34 are shown to be odd w.r.t. ^A(k) — t- — ^^(k). Thus in a particle-hole symmetric 
theory, these terms will vanish upon integration over ^A(k) and are labeled O(pha) which stands for "particle-hole 
asymmetric" . It is convenient to rewrite these matrix elements in terms of the functions 



AA(k) = [pi+^'{k)-qi-^'{k)] [yi+\k)-e[+\k) 
$A(k) = q[+^^yx{k)+p[+^'U^)ex{k) 

ei-^)(k) 



,?i+)^(k)eA(k) 



r,i-)^(k)yA(k) 



4-^\k) 



(66) 
(67) 

(68) 



where the first one, AA(k) is referred to as the Tsuneto-function [231. ^ straightforward but lengthy calculation yields 



11 



iVi3 = O(pha) (69c) 



N23 = O(pha) (69f) 



(k) ^2 -772 -4A2(k) 



^33 = -^ 4A2(k) ^69^) 

7V34 = O(pha) (69i) 



^44 = -^-^^^(k), (69j) 



where 77 = Vk ■ q. Note, that ah expressions are vahd in the whole quasiclassical hmit, i.e. for q ^ fcp and huj <^ Ep. 
For small wave numbers, as required e.g. in the Raman case, the Tsuneto and related functions AA(k), $A(k) and o''^^ 
simplify considerably. The results for such a small-q expansion can be found in appendix 1. Our further considerations 
for response and transport properties require both main results of this section: The solution of the transport equation 
in quasiparticle space, given by Eq. dJS]), will be used directly in section „ to derive the specific heat capacity in NCS 
(see Table m. While for the discussion of the gauge mode (section „), the normal and superfluid density (section,) and 
the Raman response (section p) the non-equilibrium distribution functions after an inverse Bogoliubov transformation, 
given in Eq. ([63)) and Eq. ([69)) . are necessary. 

GAUGE INVARIANCE 

The gauge invariance of our theory is an important issue, which will be discussed in the following section. Therefore, 
we determine the gauge modes and insert them into the transport equations. An integration of these transport 
equations yields a continuity equation which demonstrates the gauge invariance of our theory for Huj <^ Ep and 
q ^ fcp- For this purpose, it is very instructive to rebuild the original distribution function by combining (5n^ and 
Sn^ from Eq. ^ and Eq. ([SS)) : 



The left hand side of this equation is of the same structure as the linearized Landau-Boltzmann equation of the 
normal state. In what follows, we want to discuss the right hand side of the above equation. Note that all terms 
coupling to (5A J have vanished because of particle-hole symmetry. This means that the amplitude fluctuations of the 
order parameter do not contribute to the response in a particle-hole symmetric theory. The phase fluctuations are 
also given by Eq. : 



2 4Ai 



-A 



^A,- = -^^^Aa . (71) 
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Multiplication with the pairing-interaction V^^, and summation over k' and the band-index fi yields 



6A-{k) + J2V^^^, 



2 9 

uj — n 



6A;ik') (72) 



- 2_.^kk' 2A4k') ^^^^^ ' 

where we have introduced 69 \ = 9'^ /2 — 9x. It can be shown, that the (k)-integral over 69 vanishes identically for 
all q. Using the equilibrium gap-equation [Eq. (jl2p ] we arrive at 

^|A,|4;.^kk'4|^^(k,)|V(k)-2^T^kU' 2A,(k') ^^^^>- ^^^^ 

These arc two coupled equations (for A = ±) which determine the phase fluctuations of the order parameter (gauge 
mode) . Note that in the weak coupling BCS theory, there are only two collective excitations possible: the Anderson- 
Bogoliubov and 2A mode. In NCS, there exist two gauge modes due to the band splitting, which can be connected 
with the particle number conservation law. In addition, due to existence of a triplet fraction, there could be further 



collective excitation analogous to Leggett's SBSOS [21j modes predicted for the superfluid phases of ■^He. The latter 
should be connectable with the spin conservation law in NCS. Finally, massive collective modes with frequencies below 
2A//i may exist in NCS. It can be shown, that the right hand side of Eq. (f70| vanishes upon k and A (band) summation 
when inserting the above expressions for the gauge mode. This leads us to the following continuity equation for the 
electron density: 



^ 6nx (k) - q • ^ Vk [6nx (k) + $ a (k)^^ (k)] = . (74) 

k,A k,A 

For a conserved quantity such as the particle or charge density flk = 1, e we can identify the corresponding generalized 
density and current density 

6na = ^ak<5nA(k) (75) 

k,A 

ja = ^akVk[<5nA(k) + $A(k)(56(k)] , (76) 

k,A 

obeying the continuity equation 

uj6na - ci- ja = ■ (77) 
Therefore, we have demonstrated charge conservation and gauge invariance of the theory for hui <C Ep and q ^ /cp • 

NORMAL AND SUPERFLUID DENSITY 

The normal and superfluid density are derived in the static and long-wavelength limit (w — )■ and q — 0) . In order 
to preserve gauge invariance, gradient terms of the order 0(q) are still taken into account. The parity-projected 
distribution functions are obtained from Eq. (|63p and from Eq. (j69]): 



fe+(k) = -Mmii^) (78) 

6n- (k) = - [0, (k) - Xx (k)] 6^- (k) + rjXx (k) ^^^TT ' (79) 



where we made use of the q — >■ limit with the coherence-factors <z^(k) — ?> 1, p^(k) — ^ 0, and <I>A(k) — >■ 0a (k), 
^^(k) — )• yA(k), as well as the Tsuneto-function AA(k) — )■ (f>xO<.) — yxO<-) (see appendix 1). The combined expression 
for (5n^(k) and 6n^(k) are now inserted in Eq. (j76p to derive the supercurrent density (vertex-function a^^ — e): 

it = J2^^pd^Mp) + Mpnx(p)] (80) 

pA 



e2^VpjVpj AA(p) ( - -A+ - V6ipx 

pA 



13 



Here we used the result from section □ that 6A^ (k) / A\ (k) = i5(p\ represents the phase fluctuations of the order 
parameter. These phase fluctuations ensure gauge invariance in the above expression for the supercurrent. By 
rewriting the supercurrent as product of the superfluid density and the corresponding velocity v", we can easily 
identify 

j" = en" - (81) 
v'^ = - (--A+^V6^x) ■ (82) 

771 \ C 2 / 

Therefore, the superfluid and normal fluid density tensor read 

= ^Pi^P3>'>'iP) (83) 

pA 

< = n6,j - = ^ P^VpjyA(p) • (84) 

pA 

Thus, in this static and small~q-limit we obtain a very clear picture: The Yosida-kernel y\{\i) — —df[E\{\i)\/dE\{\i) 
generates the normal fluid density and the Tsuneto-function AA(k) gives rise to the superfluid density. 

It is important to realize, that this result can be derived in the following alternative simple way from local equilib- 
rium considerations. In terms of the Fermi-Dirac distribution function on both bands f[E\(p)\ for the Bogoliubov 
quasiparticles, the supercurrent can be written in the standard quantum-mechanical form: 

Jl = ^ + E (P)/(^A (P) + P • v^) (85) 

pA 



pA 



n5,i - — > — I I Pi > V 



V dEx{v) 

pA 



This immediately implies the definition of the normal fluid density in the form 

pA 

Thus, the results obtained with our simple local equilibrium picture are in agreement with the results in Ref. Q. 

THE SPECIFIC HEAT CAPACITY 

In order to derive the specific heat capacity, we start from an expression for the entropy of a NCS, which has to be 
written in the general form 

Ta{T) = -:^5]/[i?,(p)]ln/[i?,(p)] + {l-/[i?,(p)]}ln{l-/[i?,(p)]} 

pA 

pA 

The change of the entropy as a consequence of a temperature change 5T can then be written in the form [qI 

T5o{T) - l^ii;;,(p)57^,(p), (88) 

pA 

where the quasiparticle distribution function is given by Eq. (|45p . In the static and homogenous limit, i.e. w — > 
and q — > 0, this expression simplifies considerably to 5v\{\<i) — y\i\<i)5E\i\<i). The quasiparticle energy shift for a 
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temperature change is 5Ex{k) = {Ex{k)/T - dExik)/dT)ST for each band Q. Therefore, our result for the entropy 
change reads 



TSaiT) = l^y,(p)iJ,(p) 

pA 

^ Cv{T)6T 



Exip)-T 



dT 



6T 



and one may easily identify the specific heat capacity as 



pA 



dT 



(89) 



(90) 



An alternative way to derive the specific heat capacity employs again the concept of local equilibrium: 



V 

pA 

The change of the BQP (Fermi-Dirac) distribution function with temperature has two causes: first the direct change 
T ^ T + ST and second the change of the BQP energy with temperature through the T-dependence of the energy 
gap: 



SfiExip)) = / 



'Exip) + ^^ST^ 



dT 



I 



fEx{p) 



kB[T + 5T] J ' \ kBT 

df{Ex{p) )\ (Elx dExip 
dEx{p) 



(92) 



T 



dT 



ST 



S/a(p) 

Hence we arrive to the same result for the entropy change 



TSa{T) = l.Y.y,{p)Ex{p) 

pA 

= Cv{T)ST 



Ex{p)-T 



d Exjp) 
dT 



6T 



(93) 



and the result for the specific heat capacity is confirmed. Again, like in the case of the normal and superfluid density, 
the result for the specific heat capacity can be viewed to consist of contributions from the two bands, in the sense 
that the sum over the spin projections cr = ±1 is replaced by a sum over the pseudospin variable A = ±. 

A CASE STUDY: RAMAN RESPONSE 

In the following section we will discuss in detail the electronic Raman response for T = in NCS ^8|| . An extensive 
description of the electronic Raman effect in unconventional superconductors can be found in Ref. A Raman 
experiment detects the intensity of the scattered light with frequency-shift uj = luj — LUg , where the incoming photon 
of frequency w/ is scattered on an elementary excitation and gives rise to a scattered photon with frequ ency lus and 
a momentum transfer q. The differential photon scattering cross section of this process is given by Ref. [19[ 



5^cr ujs 2o r \ 



(94) 



with the solid angle fl and the Thompson radius vq = e'^/mc^. The generalized structure function 5-y^(q, cj) is con- 
nected through the fluctuation-dissipation theorem to the imaginary part of the Raman response function x^-y(q, w): 



^^^(q,^) 



— [l + n(^)]x"^(q,w) 



(95) 
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Here, n(a;) = [ex-p{Huj / kBT) — 1] denotes the Bose distribution. After Coulomb renormalisation and in the long- 
wavelength limit (q = 0), the Raman response function is given by the imaginary part of (see also Ref. (23| ) 

Within our notation, the unscreened Raman response is given by 



Xyyiu^) = X^^^iu^) ~ ^ (0), . • (96) 



xi°H^) = |E«pMpH. (97) 

p,(T 

( R) 

where the vertex-functions Op, bp are either 1 or the corresponding momentum-dependent Raman vertex 7 = 
that describes the coupling of polarized light to the sample. The long-wavelength limit of the Tsuneto-function 
Ap(q = 0) = 4Ap0p/ (4i?p — w^) is given in appendix 1 and since we are interested in the T — Raman response it 
is possible to perform the integration on the energy variable (see e.g. Q). Note that the second term in Eq. (|96|) 
is often referred to as the screening contribution that originates from gauge invariance. Since the ASOC leads to a 
splitting of the Fermi surface, the total Raman response is given by x^°j^^ — Sa=± ^^^^ ^77 ~ X77(^±)7 in 
which the usual summation over the spin variable a is replaced by a summation over the pseudo-spin (band) index 
A. With Eq. ([TT|) the unscreened Raman response for both bands in the clean limit [I ^ ^(0) with the mean free path 
I and the coherence length £,{T — 0)] can be analytically expressed as 




Here, Np reflect the different densities of states on both bands and (. . .)fs denotes an average over the Fermi surface. 
We consider small momentum transfers (q — > 0) and neglect interband scattering processes, assuming non-resonant 
scattering. Then, the Raman tensor is approximately given by 

where e'^'^ denote the unit vectors of scattered and incident polarization light, respectively. The light polarization 

( 

selects elements of this Raman tensor, where ^j^. can be decomposed into its symmetry components and, after a 
straight forward calculation (see appendix 2), expanded into a set of basis functions on a spherical Fermi surface. Our 
results for the tetragonal group C^y are 

00 l<k/2 

= E E "fKi ^' (100a) 

00 i<(fc+l)/2 

7s?-E E 7lf cos(4Z-2)0sin2'=0, (100b) 

k=l 1 = 1 

00 i<(fc+l)/2 
(R) J^) a;^(Al _n\A. .--.^^k , 



Tb'I-Y. E lifsm{Al-2)cj^sin''^e, (100c) 
and for the cubic group O we obtain 



fe=i 1=1 



00 l<k/2 

=J2J2 ^if sin^'^ e, (101a) 

fc=o 1=0 

7^^^, =7^''^(2-3sin2 0) + ..., (101b) 

00 ;<(fc+i)/2 



7s'?i,=E E lifcos{Al-2^,in''^9, (101c) 
fc=i 1=1 
00 ;<(fc+i)/2 

7Tf = E E 7i?sin(4Z-2)0sin2'=0 (lOld) 



fc=i 1=1 
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FIG. 2: Calculated Raman spectra for a pure triplet order parameter (i.e. tp = 0) for Bi,2 polarization of the point group 
Civ in backscattering geometry (zz). The ABM (axial) state with |dk| = dosinS is displayed in blue and the polar state with 
|dk| = rfo| cos^l in green. For a comparison, also the threshold-behavior of the Raman response for the BW state (red) with 
|dkl = do is shown. 



in a backscattering-geometry experiment (zz) [3j|. In what follows, we neglect higher harmonics and thus use only 



the leading term in the expansions of 7^^'' [34 1 



In general, due to the mixing of a singlet and a triplet component to the superconducting condensate, one expects 
a two-peak structure in parity-violated NCS, reflecting both pair-breaking peaks for the linear combination [see 
Eq. (jlip ] of the singlet order parameter ip^ (extensively discussed in Ref. [5]) and the triplet order parameter dk 
(shown in Fig. [2]), respectively. The ratio p = d/ip, however, is unknown for both types of ASOCs. 

How does the Raman spectra look like for a pure triplet p-wave state? Some representative examples, see Fig. [2j 
are the Balian-Werthamer (BW) state, the Anderson-Brinkman-Morel (ABM or axial) state, and the polar state. 
The simple pseudoisotropic BW state with dk = dgk [equivalent to Eq. ([3|) for ga = 0], as well as previous work 
on triplet superconductors, restricted on a (cylindrical) 2D Fermi surface, generates the same Raman response as an 
s-wave superconductor 17]. However, in three dimensions we obtain more interesting results for the axial state with 
dk = do{kyex — kxGy) [equivalent to Eq. ([2]) for g|| = 0]. The Raman response for this axial state in Bi and B2 
polarizations for Q = dv is given by 



7riVF7o 



(R)2 



128 
-10- 



(102) 



28 



5 + 3a;2 + 3a;'' + Sa;^ 



In 



1 



with the dimensionless frequency x = Lo/2d[). An expansion for low frequencies reveals a characteristic exponent 
[x'bi 2 ^ (1^/2^0)^] J which is due to the overlap between the gap and the vertex function. Moreover, we calculate the 
Raman response for the polar state with dk = dokze^; in this case one equatorial line node crosses the Fermi surface 



and we obtain: 




(103) 



with the trivial low frequency expansion Xbi 2 ^ w/2(io. While the pair-breaking peaks for the BW and ABM state 
were both located at w = 2dQ (similar to the Big polarization in the singlet d-wave case, which is peaked at 2Ao), for 
the polar state this peak is significantly shifted to lower frequencies (w = 1.38do)- 

Let's turn to the predicted Raman spectra for the tetragonal point group G = C^y. In Fig.|3]we show the calculated 
Raman response using Eq. ^ with g|| = 0. This Rashba-type of ASOC splits the Fermi surface into two bands; while 
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FIG. 3: Calculated Raman spectra X77 (A- ) [blue] andx"7(A+) [red] for Ai (solid lines) and for Bi, 2 (dashed lines) polarizations 
for the point group C4„. We obtain the same spectra for the Bi and B2 symmetry. The polar diagrams in the insets demonstrate 
the four qualitative different cases for the unknown ratio p = 



on the one band the gap function is Ak = V' (1 + p|7kl) = A+, it is A_ = ?/) (1 — p|7k|) on the other band. Thus, 
depending on the ratio p = four different cases (see polar diagrams in the insets) have to be considered: (a) no 
nodes; (b) one (equatorial) line node (A_ band); (c) two line nodes (A_ band); and (d) two point nodes on both 
bands. Since the Raman intensity in NCS is proportional to the imaginary part of 

x;°*^' = X77(A-) + X77(A+), (104) 

it is interesting to display both contributions separately (blue and red lines, respectively). Although (except for V-" = 0) 
we always find two pair-breaking peaks at 

^ = |1±P| (105) 

we stress that our results for NCS are not just a superposition of a singlet and a triplet spectra. This is clearly 
demonstrated in Fig. [3Ja), for example, in which we show the results for a small triplet contribution (p = 1/2). For 
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FIG. 4: Calculated Raman spectra x-yyi^-) [blue] and X'n{^+) [red] for E (solid lines), T2 (dashed lines) and Ai [dotted line, 
only in (d)] polarizations for the point group O. The insets display the point and line nodes of the gap function A_. 



x"^(A_) we find a threshold behavior with an adjacent maxinium value of x'^^^(A_) = -^p 7o ^ "^"^ l^\/v ^ 1- 
In contrast for x"').(A+) a zero Raman signal to twice the singlet contribution followed by a smooth increase and 
a singularity is obtained [SS*]. In the special case, in which the singlet contribution equals the triplet one (p — 1), 
the gap function A_ displays an equatorial line node without sign change. This is displayed in Fig. Because of 
this nodal structure and strong weight from the vertex function (oc sin^ &) , many low energy quasiparticles can be 
excited, which leads to this square-root-like increase in the Raman intensity. In this special case the pair-breaking 
peak is located very close to elastic scattering (w = 0.247/1). In Fig. [Sjc) the gap function A_ displays two circular 
line nodes. The corresponding Raman response for p > 1 shows two singularities with different low frequency power 
laws [xsi Jy^-) ^ and x'bx 2^^+) (w/2-!/' — 1)^^/^]. Finally, for p 1 one recovers the pure triplet cases (d) 

which is given analytically by Eq. (I102[) . 

The Raman response for the point group O, using Eq. ([3]), is shown in Fig. |4l We again consider four different 
cases: (a) no nodes; (b) six point nodes (A_ band); (c) six connected line nodes (A_ band); and (d) 8 point nodes 
(both bands) as illustrated in the insets. Obviously, the pronounced angular dependence of |7j(.| leads to a strong 
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polarization-dependence. Thus we get different peak positions for the E and T2 polarizations in x"^(A+). As a 
further consequence, the Raman spectra reveals up to two kinks on each band (+,— ) at 

^ = |l±p/4| (106) 

and 

^^\1±P\. (107) 

Interestingly, the T2 symmetry displays only a change in slope at uj/2'4> — |1 +p| instead of a kink. Furthermore, no 
singularities are present. Nevertheless, the main feature, namely the two-peak structure, is still present and one can 
directly deduce the value of p from the peak and kink positions. Finally, for 1 one recovers the pure triplet case 
(d), in which the unscreened Raman response is given by 



" (fl)2 |7kP 

^ \ " v/(^/2d+|7kl)(W2d-|7kl)/ PS 



Clearly, only the area on the Fermi surface with ijj/2d > \ contributes to the Raman intensity. Since \ E [0, 1] 
has a saddle point at = 1/4, we find kinks at characteristic frequencies uj/2d =1/4 and uj/2d = 1. In contrast to 
the Rashba-type ASOC, we find a characteristic low energy expansion cx {uj/2d)'^ for both the Ai and E symmetry, 
while oc {ijj/2d)^ for the T2 symmetry. Assuming weak coupling BCS theory, we expect the pair-breaking peaks (as 
shown in Fig. S]) for Li2Pda;Pt3_:rB roughly in the range 4cm~^ to 30cm~^. 



CONCLUSION 



In this chapter, we derived response and transport functions for noncentrosymmetric superconductors from a kinetic 
theory with particular emphasis on the Raman response. We started from the generalized von Neumann equation 
which describes the evolution of the momentum distribution function in time and space and derived a linearized 
matrix-kinetic (Boltzmann) equation in w-q-space. This kinetic equation is a 4 x 4 matrix equation in both particle- 
hole (Nambu) and spin space. We explored the Nambu-structure and solved the kinetic equation quite generally 
by first performing an SU(2) rotation into the band-basis and second applying a Bogoliubov-transformation into 
quasiparticle space. Our theory is particle-hole symmetric, applies to any kind of antisymmetric spin-orbit coupling, 
and holds for arbitrary quasiclassical frequency and momentum with hio ^ E-p and |q| <C k-p. Furthermore, assuming 
a separable ansatz in the pairing interaction, we demonstrated gauge invariance and charge conservation for our theory. 
Within this framework, we derived expressions for the normal and superfiuid density and compared the results in the 
static and long-wavelength limit with those from a local equilibrium analysis. The same investigations were done for 
the specific heat capacity. In both cases we recover the same results, which validates our theory. 

Finally, we presented analytic and numeric results for the electronic (pair-breaking) Raman response in non- 
centrosymmetric superconductors for zero temperature. Therefore we analyzed the two most interesting classes of 
tetragonal and cubic symmetry, applying for example to CePtaSi {G = C4v) and Li2Pd2,Pt3_a;B {Q = O). Accounting 
for the antisymmetric spin-orbit coupling, we provide various analytic results such as the Raman vertices for both 
point groups, the Raman response for several pure triplet states, and power laws and kink positions for mixed-parity 
states. Our numeric results cover all relevant cases from weak to strong triplet-singlet ratio and demonstrate a charac- 
teristic two-peak structure for Raman spectra of non-centrosymmetric superconductors. Our theoretical predictions 
can be used to analyze the underlying condensate in parity-violated noncentrosymmetric superconductors and allow 
the determination of the unknown triplet-singlet ratio. 

We thank M. Sigrist for helpful discussions. 



Appendix 1: Small q— expansion 

For small wave numbers, i.e. q — > 0, the Tsuneto and related functions, which play an important role in the 
matrix-elements Nij [see Eq. (|69p ]. will simplify considerably. Taking into account terms to the order 0(?7^) with 
77k = Vk • q, we obtain the well-known expression for the Tsuneto-function [l6| 



hm A,(k) = -4Al(k) . ^ ^ ^ , \ ' (109) 
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where 



is the derivative of the electron distribution function in the band A and 

5/[i?A(k)] 1 1 



yA(k) = 



ai?A(k) 4A:BT,osh2(§M) 



(110) 



(111) 



is the derivative of the quasiparticle distribution function. 

The following limits are also of interest: the homogenous limit (q = 0), e.g. for the Raman response and the static 
limit (w = 0), used in local equilibrium situation 



AA(k,q = 0) = 



4A!(k)gA(k) 



lim lim AA(k) 

w^Oq-S-O 



</.A(k) - 2/A(k) . 

For the following small q-expansion we omitted the band-label A for better readability: 
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where /^"^ denotes the nth derivative of /(-Ek) with respect to E'k- Furthermore we find the following 
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expansions: 
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The ten products of coherence-factors in Eq. (j64l) have the following explicit form: 
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and the small q-limit of each coherence-factors reads: 
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Appendix 2: Derivation of the Raman vertices 

In order to derive the relevant expressions for the polarization-dependent Raman vertices, we start from a general 
dispersion relation for tetragonal symmetry (Ciy) 



Ck = a^^^ [cos(nfca;a) + cos(nfcj,a)] cos(rfc2c) 

n—l r— 

oo oo 

+ ^ cos{nkxa) cos{nkya) cos{rkzc) 

n=0 r=0 

oo oo oo 

+ c'^^rn,r [cos{nkxa) cos{mkya) + cos{mkxa) cos(rifcj^a)] cos(rfc2c) 



(118) 



n—l m— 1 r— 
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and for the cubic symmetry (O) 

oo 

eic = [cos{nkxa) + cos{nkya) + cos{nkzC)] 

n=l 

oo 

+ cos{nkxa) cos{nkya) cos{rkza) 

[cos{mkxa) cos{mkya) cos{nkza) 



(119) 



n=0 



oo n — 1 
n—1 m— 1 



+ cos{mkxa) cos{nkya) cos{mkza) + cos{nkxa) cos{mkya) cos(mfcza)] 

oo n—1 m— 1 

+ ^ dS,m,r [cos{nkxa) cos(mfcyo) cos(rfc2a) 

n— 2 m— 1 r— 

+ cos{nkxa) cos{rkya) cos{mkza) + cos{mkxa) cos{nkya) cos{rkza) 
+ cos{rkxa) cos{nkya) cos{mkza) + cos{mkxa) cos{rkya) cos{nkza) 
+ cos{rkxa) cos{mkya) cos{nkza)] . 

Time reversal symmetry allows only for even functions of momentum k in the energy dispersion. Furthermore the 

dispersion must be invariant under all symmetry elements of the point group Q of the crystal. For small momentum 
transfers and nonresonant scattering, the Raman tensor is given by the effective-mass approximation 



7(k) =m^ 



* h?dkidkj ^ 



(120) 



where e'^'^ denote the unit vectors of the scattered and incident polarization light, respectively. 
The light polarization vectors select elements of the Raman tensor according to 



e -Tic 



(121) 



where the Raman tensor can be decomposed into its symmetry components and later expanded into Fermi surface 
harmonics: 
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(123) 



Here we have omitted all non-Raman active symmetries such as A2g. The vertices A^^'' and A^^"* are equal up to 
some constants determined by the band structure, and the vertices for E^^^ and -E^^' in dy differ only by a rotation 
of the azimuthal angle (p by 7r/2. Since this rotation is an element of the corresponding point groups, these vertices 



are identical, too. The same holds for T^^\ T^"^^ and T.^-'> . Therefore the upper indices will be omitted in the 
following (whenever possible). For the tetragonal group €4^ the Ai, Bi, B2 and E symmetries are Raman active in 
backscattering geometry. Relevant polarizations for this group are: 
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The cubic group O reveals three Raman active symmetries, namely Ai, (E^^^ E^^^), and T2 (still assuming backscat- 
tering geometry). The relevant polarizations are: 
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Here, we have defined the unit polarization vectors x' = (x + y)/-\/2 and y' = (x — y)/-\/2- L and R denote left and 
right circularly polarized light with positive and negative helicity, respectively (e^ — (x + iy)/-\/2, = (x — iy)/\/2)- 
Note that in a backscattering configuration the polarization vectors e^''^ are pinned to the coordinate system of the 
crystal axes. Therefore some caution is advised when choosing the proper helicity for the scattered polarization vector 
e'^. Although the Raman vertices E'-'^' and E^^^ seem to look completely different, the Raman response turns out to 
be exactly the same. From a tight-binding analysis we obtain the same (band-structure) prefactors for both vertices, 
thus 7^ and \/37jf generate both the same Raman response. Note that it is not possible to measure Ai and E^^) 
independently in backscattering geometry with the crystal c-axis aligned parallel to the laser beam. 

The Raman vertices are extracted from the band structure by comparing the symmetry components of the Raman 
tensor with the second derivative of the energy dispersion. This can be done by solving a set of 6 coupled linear 
equations - the 6 equations correspond exactly to the 6 free components of the symmetric tensor of inverse effective- 
mass and to the 6 symmetry elements (vertices) to be determined. Finally we make a series expansion in k, in order 
to get the angular dependence of the vertices on the Fermi surface. Our results for the tetragonal point group C4„ 
are 

00 l<k/2 

IM-Y. E 7S^cos4Z0sin^'=0 (126a) 

fc=0 (=0 
00 i<(fc+l)/2 

7i?=E E 7i?cos(4/-2)0sin2'=0 (126b) 
fc=i 1=1 

00 i<(fc+l)/2 

7i?-E E 7i?sin(4/-2)0sin2'=0 (126c) 



00 00 



' E E ^kf sin(2/ - 1)(^ sin 2k9 (126d) 

k=l 1 = 1 



and for the cubic point group O we obtain 

00 l<k/2 

=Y.J2 ^if 4/0 sin^'^ e (127a) 



fc=o ;=o 



7),^', =7r(2-3sin^0) + ... (127b) 

00 ;<(fe+i)/2 

^^^l = E E ^if - 2)0 sin^'^ e (127c) 

k=l 1 = 1 

00 ;<(fe+i)/2 

7Tf = E E ^if «i"(4/ - 2)0 sin^'^ (127d) 

k=l 1=1 

in a backscattering-geometry experiment (z'z) . 
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